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Using 53 pb“^ of e'’"e“ annihilation data taken at y/s = 3.686 GeV, a comprehensive study has 
been made of the radiative decays of samples of 5.1 million J/%p and 24.5 million ’tp{2S) into pairs 
of pseudoscalar mesons, 7r''"7r“, K~ , K%K%, and rjrj. Product branching fractions for the 

radiative decays of J/'ijj and ip{2S) to scalar resonances /o(1370,1500,1710, 2100, and 2200), and 
tensor resonances /2(1270,1525, and 2230) have been determined, and are discussed in relation to 
predicted glueballs. For tp{2S) radiative decays the search for glueballs has been extended to masses 
between 2.5 GeV and 3.3 GeV. 


I. INTRODUCTION 

Quantum Chromodynamics (QCD) is presently the ac¬ 
cepted theory of strong interactions. Its predictions for 
hadronic phenomena at high energies have received re¬ 
markable experimental confirmation, and it is an inte¬ 
gral part of the Standard Model of particle physics. In 
contrast, at lower energies the success of QCD has been 
comparatively limited. This is mainly due to the fact 
that while at high energies reliable predictions can be 
made in perturbative QCD (pQCD) models, at low ener¬ 
gies pQCD predictions are unreliable, and lattice calcu¬ 
lations with the required precision are in many cases not 
yet available. 

In 1972, Fritzsch and Gell-Mann [I| pointed out that 
“if the quark-gluon field theory indeed yields a cor¬ 
rect description of strong interactions, there must exist 
glue states in the hadron spectrum”, and Fritzsch and 
Minkowski Q presented a detailed discussion of the phe¬ 
nomenology of the spectrum of what they called “glue- 
states”. The “glue-states”, which we here call “glue¬ 
balls” , are bosons with baryon number and isospin equal 
to zero, and they are SU (3) flavor singlets “to the extent 
to which SU{3) breaking effect can be neglected” 0. 

Ever since 1985, when ?7(1440) was claimed as the first 
glueball discovered Q, numerous experimental searches 
for glueballs have been reported, and numerous claims 
and counterclaims abound. Many reviews of QCD ex¬ 
otics, including glueballs, have been published. Two of 
the latest comprehensive reviews are by Klempt and Za¬ 
itsev 0 and Ochs 

Fritzsch and Minkowski Q have detailed color-singlet 
states of two and three gluons allowed by conservation 
rules, which are identical to those for states of two and 
three photons. The two-gluon states have charge con¬ 
jugation, C = -1-1, while the three gluon states can 
have C = -1-1 or —1. The lowest mass states have S- 
waves between the gluons, and are two-gluon states with 
= 0++, 2++, and 0 and three-gluon states with 
= 1++, l+“, 1 , and 3 . 

Predictions for glueball spectra have been made in 
bag models QCD-based potential models Q, QCD 


sum rules 0, and lattice gauge calculations dSl. Lat¬ 
tice gauge calculations have supplanted most other pre¬ 
dictions. The lattice calculations of Refs. 10,11, and 
12 were done in the quenched approximation with no 
quarks, and their predictions are in general agreement 
with each other. The first unquenched lattice calcula¬ 
tions with 2-1-1 flavor sea-quarks have been reported by 
the UKQCD Collaboration [l^. They obtain results for 
the masses of the 0++, 0 and 2++ glueballs which 
are in good agreement with the lattice predictions of 
Refs. 10 and 12 in the quenched approximation, and con¬ 
clude that they “find no evidence of strong unquench¬ 
ing effects”. Thus, the results of Ref. 12, reproduced 
in Fig. [U can be considered as representative of lattice 
predictions so far. The order of the lowest lying glue¬ 
balls is as anticipated by Fritzsch and Minkowski Q, 
with the prediction for the first scalar (O^+j and tensor 
(2++) glueball masses are G(0j''^) = 1710(50)(80) MeV, 
= 2390(30)(120) MeV. These are followed by five 
heavier glueballs with masses less than M{ip{2S)). These 
are 0-+(2560(35)(120) MeV), 1+-(2980(30)(140) MeV), 
2-+(3040(40)(150) MeV), 3+-(3600(40)(170) MeV), and 
3++(3670(50)(150) MeV). In the present paper, our glue¬ 
ball search is limited to masses less than 3.3 GeV, below 
the masses of the Xcj resonances of charmonium. 

Early searches for glueballs were misled by the wish¬ 
ful expectation of “pure” glueballs. With the realization 
that the predicted glueballs have masses close to those of 
normal qq mesons with the same quantum numbers with 
which they must mix, the more realistic searches aspire to 
identify states which contain substantial glueball content 
in their wave functions. Searches for such states have 
been based mainly on three expectations. These are: (a) 
for a given , the glueball state is “supernumary”, i.e. 
a state in excess of the expected number of qq states, 
e.g., the nonet of 0++ SU{3) states, (b) a pure glueball 
state does not decay preferentially to states containing 
up, down, or strange quarks, i.e., “pure” glueball decays 
to pseudoscalar pairs should have equal branching frac¬ 
tions for each available charged state. Thus decays to tttt, 
KK, rjr], rj'rj' and r]r]', should have branching fractions 
proportional to 3 : 4 : 1 : 1 : 0, and (c) certain decays 
are particularly favorable for populating glueballs, e.g., 
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pp annihilation, and radiative decays of vector states of 
quarkonia, J/'ip, to jgg ..where g stands for glu¬ 
ons. Numerous experimental attempts have been made 
to exploit these expectations, but they have been gener¬ 
ally limited in their scope. 

An extensive review of searches for glueballs in pp anni¬ 
hilations, central production, and hadronic and radiative 
decays of vector states of charmonium and bottomonium 
has been made by Klempt and Zaitsev Q. Radiative 
decays of have been studied by Crystal Ball 0 , 
Mark III [H [Hi, DM2 0 , and BES 0^13 ■ Among 
the latest of these are BES measurements of J/p) —>■ 
7 ( 7 r+ 7 r“, 7 r° 7 r°) [H| and J/'ip —>■ 'y{K^K~,KgKg) [H 
with 58 million directly produced J/ij: in e+e“ —>■ J/ip. 
The invariant mass spectra in the decays of these directly 
produced contain large backgrounds due to ISR pro¬ 
duction of 7 r^ 7 r“ 7 r'^ and KKtt^ with one missing photon 
which simulate 771 “’'tt” and "/K'^K~ decays. More re¬ 
cently, BES III has reported a study of J/ijj ^ 77777 with 
2.28 X 10® directly produced J/t/; 0- No studies of ra¬ 
diative decays of 'ip{2S) have been published, and only 
an unpublished report with 14 million il’{2S) by BES ex¬ 
ists 1 ^ . 

In the present paper we present results of a study of the 
radiative decays of the charmonium resonances, 5.1 mil¬ 
lion J/ 7 / 1 , and 24.5 million p^{2S) resonances populated 
in e“'"e“ annihilation, and decaying into pairs of pseu¬ 
doscalar, TTTT, KK, and 7777 . We determine the masses, 
widths, and product branching fractions of the reso¬ 
nances observed, and discuss their possible relations to 
glueballs. Since the J/ijj used in the present investigation 
are produced by the exclusive decay ip{2S) —?> tt+tt” J/t/j, 
the invariant mass spectra we observe from the J/ip 
radiative decays are free from the ISR-produced back¬ 
grounds. The invariant mass spectra for ip{2S) decays 
contain, of course, background contribution from ISR. 

The paper is organized as follows. In Sec. II we de¬ 
scribe the event selection criteria. In Sec. Ill we describe 
Monte Carlo event simulations and use them to deter¬ 
mine instrumental resolutions and efficiencies. In Sec. IV 
we present the Dalitz plots for J/t/j and '0(25') radiative 
decays to pseudoscalar pairs. In Sec. V.A we present the 
TTTT and KK invariant mass distributions for the radia¬ 
tive decays of T/ 0 , and discuss the resonances observed 
in them. In Sec. V.B we present the similar mass distri¬ 
butions and discussion for the radiative decays of tjj{2S). 
In Sec. V.C we present the results of our search for the 
putative tensor glueball ^(2230). In Sec. VI we summa¬ 
rize and compare our results for J/pj and 0(25') radiative 
decays to pseudoscalar pairs for the observed resonances 
with masses below 2.5 GeV. In Sec. VILA we present the 
invariant mass spectra for 0(25) radiative decays to tttt, 
KK, and 7777 decays in search for potential glueball can¬ 
didates in the mass range 2.5 — 3.3 GeV, and establish 
upper limits for any resonance enhancements. We ob¬ 
serve strong decays of Xco( 0 ’*’''’) and Xc 2 ( 2 ’^'*') resonances 
of charmonium in our data for 0(25) radiative decays. 
Because the decays of these resonances are necessarily 
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FIG. 1. Quenched lattice predictions of glueball spectrum by 
Chen et al. [13] ■ The states are grouped by their parity and 
charge conjugation indicated in the abscissa. 


TABLE I. Summary of the requirements for the number of 
reconstructed particles for the reactions investigated in this 
paper. 


Decay mode 

0 tracks 

0 showers 

Other 

0(25) —>■ 777 ■^77“ 

2 

> 1 


-iK+K- 

2 

> 1 


0 

C 

0 

> 5 

2 77° 

77777 

0 

> 5 

277 

'yKiKl 

> 4 

> 1 

2Kl 

J/tp 777+77 

4 

> 1 


-fK+K- 

4 

> 1 


0 

0 

2 

> 5 

2 77° 

77777 

2 

> 5 

277 

'jKlKl, 

> 6 

> 1 

2 Ag 


mediated by two gluons, just as the decays of two-gluon 
glueballs, it is instructive to examine and compare their 
characteristics. We do so in Sec. VII.B in terms of the 
branching fractions measured using the same data as we 
use in the present paper, and published as Ref. [13. In 
Sec. VIII we discuss systematic uncertainties in our re¬ 
sults. In Sec. IX we summarize the results of the present 
investigation for the decays of known and proposed reso¬ 
nances into pseudoscalar pairs, and their implications for 
scalar and tensor glueball candidates. 









3 




o 

O 




FIG. 2 . Two-photon invariant mass distributions for radia¬ 
tive decays of (top panel) J/^ and (bottom panel) ip{ 2 S), 
with at least 5 photons in the final state. In both panels, the 
top histograms show invariant mass distribution for all two- 
photon combinations in the selected events. The middle his¬ 
tograms show distributions for all two-photon combinations 
in the events that pass the 4 C energy-momentum kinematic 
fit with xlc < 25 . The shaded solid histograms show the dis¬ 
tributions of events which contain two 7 r° or two rj candidates 
which best fit and T] masses. Signals from 7r° 77 and 

77 —>■ 77 decays are clearly visible in all three histograms. 


FIG. 3 . Two-pion invariant mass distributions for radiative 
decays of (top panel) J/ip and (bottom panel) 'tp{ 2 S) with 
at least two Kg Tr~^ 7 r~ candidates reconstructed. In both 
panels the top histograms show all Kg —>■ candidates 

in the selected events. The dashed histograms show distribu¬ 
tions for all Kg —> 7 r~^Tr~ candidates in the events that pass 
the 4 G energy-momentum kinematic fit with xic < 25 . The 
shaded histograms show the distribution for the best two Kg 
candidates in the event, as determined by the kinematic fit. 
The mass peaks corresponding to the selected Kg candidates 
are clearly visible in all three histograms. 


II. EVENT SELECTIONS 

We use 53 pb“^ of e+e“ annihilation data taken at 
^/s = 3.686 GeV at the CESR e+e“ collider at Cor¬ 
nell University. It corresponds to 24.5 million f/’(25') 
produced, and 5.1 million J/ij} tagged by the decay 
ip{2S) —>■ 7 r+ 7 r“ J/r/j. We also use 21 pb“^ of e+e“ an¬ 
nihilation data taken at ^/s = 3.672 GeV to study back¬ 
grounds from off-resonance decays. The decay products 
were analyzed by the CLEO-c detector, which has been 
described elsewhere in detail [ 2 ^. Briefly, it consists of 
a CsI(Tl) electromagnetic calorimeter, an inner vertex 
drift chamber, a central drift chamber, and a ring imag¬ 
ing Cherenkov (RICH) detector, all inside a supercon¬ 
ducting solenoid magnet providing a nominal 1.0 Tesla 
magnetic field. The acceptance for charged and neu¬ 
tral particles is |cos0| < 0.93. Charged particle reso¬ 
lution is CTp/p = 0.6% @ 1 GeV/c. Photon resolution 


is (Te/E = 2.2% @ 1 GeV, and 5% @ 100 MeV. The 
detector response was studied using a GEANT-based 
Monte Carlo (MG) simulation including radiation cor¬ 
rections [^ . 

The particle multiplicity requirements for the vari¬ 
ous decay modes are given in Table ID Charged parti¬ 
cle tracks were accepted in the full range of the drift 
chamber’s acceptance, | cos0| < 0.93 (0 is the polar angle 
with respect to the positron beam), and were selected by 
standard CLEG quality requirements. The total charge 
in the event was required to be zero. Photon candi¬ 
dates were accepted in the “Good Barrel” and “Good 
Endcap regions” of the calorimeter (|cos0| < 0.81 and 
I cos0| = 0.85 — 0.93, respectively), and were required to 
contain none of the few noisy calorimeter cells, and to 
have a transverse energy deposition that was consistent 
with that of an electromagnetic shower. 

The 7 r° and 77 were reconstructed in their decay to two 
photons. Pairs of photon candidates were kinematically 
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Recoil Mass against n-^n (MeV) 



FIG. 4 . Mass distribution of events recoiling against 7r"''7r“ in 
tp{ 2 S) —>■ J/ip decay (a) before and (b) after 4 C kine¬ 

matic fit. The fitted mass of the J/ip peak is found to be 
M( J/V') = 3097.13 ± O.Ol(stat) MeV. 


FIG. 5 . Distributions of |cos0(/o,2)| for / 2 ( 1270 ) —>■ 7r’''7r“ 
and /o(1710) —>■ K^K~ decays. The point show the distribu¬ 
tions from data, and the shaded histograms show the results 
of MG simulations. 


fitted to the known and rj masses, and the pairs were 
accepted if they had a mass < 3cr of the particle mass. 
Invariant mass distributions at different stages of these 
selections are shown in Fig. [2l 

The Kg were reconstructed in their decay to two 
charged pions. Pairs of oppositely charged pions were 
kinematically fitted to originate from a common vertex, 
and this common vertex was required to be displaced by 
> 3(7 from the e“'"e“ interaction point. Invariant mass 
distributions at different stages of these selections are 
shown in Fig. O 

Charged pions and kaons were identified using infor¬ 
mation from both dE/dx and RICH measurements. For 
the pion and kaon tracks with momenta > 600 MeV, the 
variable was constructed, combining information 

from dE/dx measurements, and the RICH log- 

likelihood, , for the particle hypotheses i and j 




■ ^dE ldx\2 




RICH 


tRICH 

^3 


To identify pions, it was required that ALt^^k < 0. For 
kaons, it was required that AL,r,if > 0. For pion and 
kaon tracks with momenta < 600 MeV, for which RICH 
information was not available, it was required that the 
dE/dx of the track be consistent within 3a of the re¬ 
spective mass hypothesis, and that it be more consistent 
with its mass hypothesis than with the alternate mass 


hypothesis. 

For the 77r“'"7r“ final state, electron contamination was 
removed by rejecting events with Ecc/p > 0.85, and 
muon contamination was removed by rejecting events 
which have tracks matched to any hits in the muon cham¬ 
bers. 

Since we study resonances with masses < 2.8 GeV for 
J/4> decays and < 3.3 GeV for t/>{2S) decays, we safely 
require E^ > 80 MeV for the “radiative” photon candi¬ 
date. 

We obtain our J/tp sample from tjj{2S) by tagging 
events with two oppositely charged tracks, which were 
assumed to be pions, to have a recoil mass in the range 
M^J/pi) ± 10 MeV. The mass distribution of events re¬ 
coiling against tt+tt” in the decay ip{2S) —^ J/ip is 

shown in Fig.|4l[a) before 4C kinematic fit, and Fig. 4(b) 
after 4C fit. The strong peak corresponds to J/ip events. 

To reconstruct the full ip{2S) and J/ip events, we make 
a 4C energy/momentum kinematic fit, constraining the 
sum of the 4-vectors of all final state particles to be equal 
to the 4-vector of the e+e" interaction. For J/ip decays, 
the J/ip decay products were additionally constrained to 
have M{J/ip). In all cases the final fit was required to 
have xic <25. When constructing the invariant masses 
from the radiative decays, the 4-vectors resulting from 
these fits were used. 
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FIG. 6. Dalitz plots for (top two panels) J/ip —>■ 777and 777^77°, and (bottom two panels) J/t/i —>■ 'yK'^K and ')K%K%. 
The dashed lines indicate the regions used to reject J/t/i —>■ pir and K*K decays in the 777’^77“ and ■yK^K~ decays, respectively. 


III. MONTE CARLO SAMPLES AND 
EFFICIENCIES 

Monte Carlo samples of simulated events were gener¬ 
ated using EvtGen [ 2 ^, and then passed through the 
standard detector simulation and reconstruction rou¬ 
tines, to determine mass resolutions and experimental 
efficiencies. Monte Carlo samples of 200k events were 
generated for each pseudoscalar pair decay chain for the 
light quark resonances seen in the data and the presumed 
glueball candidates in the mass region 2500 — 3500 MeV. 
To investigate background contributions from other de¬ 
cays, such as J/ 7/7 —>■ pTr,K*K, Monte Carlo samples of 
5 million events were generated. When considering radia¬ 
tive decays to a particular resonance (e.g. ip —>■ 7 / 0 , 2 ), 
the expected 1 + a cos^ 0 angular distributions for the 
corresponding El transitions were used, with a = 1 for 
decays to a 0++ state and a = 1/13 for decays to a 
2^'^ state. The angular distributions for scalar and ten¬ 
sor mesons decaying into two pseudoscalars were simu¬ 
lated using the helicity amplitude formalism of EvtGen. 
Fig. [S] illustrates excellent agreement between the data 


and simulated angular distribution for the well-separated 
resonances /2(1270) and /o(1710). 

The efficiencies obtained for decays to individual reso¬ 
nances are shown in Table HIl For J/tp decays to charged 
final states the efficiencies range from 15% to 20%, and 
for neutral final states they range form ~ 13% to 16%. 
Corresponding efficiencies for ip{2S) decays are approxi¬ 
mately a factor two larger. 

Monte Carlo simulations give instrumental Gaussian 
resolution widths for masses < 2500 MeV as cr = 3 — 
4 MeV, for tt+tt”, K^K~, and KgKg, 9 — 11 MeV for 
77 ^ 77 °, and 7 — 8 MeV for 7777 . 


IV. DALITZ PLOTS FOR J/ip AND ip{ 2 S) 
RADIATIVE DECAYS 

J/xjj ^ 77777: 

Dalitz plots for the radiative decays of J/tp —>■ yTr+Tr” 
and 77 ^ 77 ° are shown in the top two panels of Fig. [51 In the 
777 + 77 “ plot there are horizontal and vertical bands corre¬ 
sponding to p(770). As confirmed by Monte Garlo simu- 
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TABLE II. MC-determined efficiencies for the reactions 
V'(2S) 7 i?, PS, PS, and V'(2S) ->■ -k+tt-J/iP, J/i/> -)■ 

yR, R —^ PS, PS, for selected intermediate resonances R. 


Resonance 

Efficiency, (%) 

X 

J/ip ^ yX 

ip{2S) yX 

X —>■ 

/2(1270) 

19.8 

36.4 

/o(1500) 

18.5 

35.9 

/o(1710) 

17.1 

34.5 

/o(2100) 

15.0 

34.1 

2500 - 3300 

— 

37.5 - 29.4 

X —>■ 7r“ 

77 


/2(1270) 

13.1 

24.5 

/o(1500) 

13.2 

25.3 

/o(1710) 

13.8 

25.7 

/o(2100) 

14.8 

26.9 

2500 - 3300 

— 

27.4 - 30.0 

X-^ K+K- 

/o(1370) 

/2(1525) 

23.3 

36.1 

23.2 

35.2 

/o(1710) 

22.4 

35.6 

/o(2200) 

18.6 

36.5 

2500 - 3300 

— 

40.6 - 40.8 

K'/ 

—>■ 7r'*’7r~ 


/o(1370) 

/2(1525) 

14.7 

23.1 

15.0 

23.0 

/o(1710) 

15.0 

23.6 

/o(2200) 

15.9 

24.1 

2500 - 3300 

— 

27.3 - 27.5 

X —>■ 777, 7 —>■ 77 

/o(1710) 

11.3 

18.2 

2500 - 3300 

— 

17.6-16.0 


lations, these arise from decay of J/tp with one pho¬ 
ton from 7r° missing. In subsequent analysis we remove 
these by rejecting events with M‘^{Tr^y) < 1.0 GeV^ . 
We do not remove tt+tt” events from decays which 
populate the diagonal band in the plot in order to be 
able to take the effect of their tails in subsequent fits to 
the 7r+7r“ mass spectrum. Fig. [7{a) shows that the tail 
of the residual p peak does not contribute significantly to 
the 7r+7r“ mass spectrum in the M(7r+7r“) > 1.2 GeV. 

As expected, the p(770) contributions are absent in the 
7r°7r° Dalitz plot. 

The only other feature which is identifiable in these 
Dalitz plots is the excitation of /2(1270) as a diago¬ 
nal band with J = 2 characteristic angular distribu¬ 
tion which produces enhancements near the maxima of 
M^{TTy). 

J/ip jKK: 

Dalitz plots for the radiative decays of J/ip — >■ yK^K~ 
and KgKg are shown in the middle two panels of Fig. |6l 
The bands due to p(770) excitation seen in the 77r+7r“ 
plot are , of course, absent in these plots. Instead, 
there is evidence for contributions from the decay J /ip ^ 
iF+iF“*(892) —7> iF+iF“7r°(-|-c.c.), with one photon from 
TT^ missing. We remove these events in the horizontal 



FIG. 7. (a) Distributions of 7r'*’7r~ invariant mass for J /ip ^ 
77r"*’7r“ decays for data (solid histogram) and J/ip pn 
MC decays (shaded histogram), illustrating the small residual 
non-p background from pn decay, (b) Distributions of K~ 
invariant mass for J /ip ^ yK^ K~ decays for data (solid his¬ 
togram) and J /ip KKtx^ MC decays (shaded histogram). 


and vertical bands by rejecting events with M‘^{K^y) < 
1.0 GeV^. Monte Carlo simulation, shown in Fig. [3(b), 
confirms that the residual background, which is too small 
to be discernible, is completely negligable. Excitation of 
72 ( 1525 ) with its characteristic J = 2 angular distribu¬ 
tion is identifiable as a diagonal band, followed by the 
nearly uniform population of /o(1710). A weakly pop¬ 
ulated diagonal band corresponding to /o(2200) is also 
discernible. 

The J/ip ^ yKgKg Dalitz plot is similar to the of 
J /ip ^ yK^K~ plot but has smaller yield. 

ip{2S) —)• 77r7r, 'jKK: 

The four Dalitz plots for the tttt and KK decays of ip{2S) 
are shown in Fig. [51 As is well known, all hadronic and 
radiative decays of ip{2S) are suppressed by a factor ^ 7 
(the so-called 13% rule) compared to those of J/ip. As 
a result all plots in Fig. [8] have smaller yields than in 
the corresponding plots in Fig. [51 . Because of the well 
known anomalously small branching fractions for the de¬ 
cays ip{2S) —>■ pTT and K*K, the p bands and the K* 
bands seen in J/ip decays are not seen in the correspond¬ 
ing ip{2S) Dalitz plots. However, there is a strong diag- 
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FIG. 8. Dalitz plots for (top two panels) —>■ yTr'^yr and yvr'^Tr®, and (bottom two panels) i/’( 2 <S') —> jK'^K and 'yK%K%. 


onal band for p(770) in the '0(25') —j> 77 r+ 7 r“ plot and 
a narrow band corresponding to 0 ( 1020 ) in 
These arise from ISR contributions as we show later. 

The strong diagonal bands due to the population of xo 
and X 2 states of charmonium are visible in all four Dalitz 
plots for 0(25) radiative decays. 


V. INVARIANT MASS DISTRIBUTIONS FOR 
J /0 AND 0 ( 25 ) RADIATIVE DECAYS 

In order to draw quantitative conclusions from our 
measurements, we make projections from the Dalitz plots 
for the invariant masses for the tttt, KK, and 7777 fi¬ 
nal states. An overall presentation of these is made in 
Fig. [3 Several gross features are apparent in this figure. 
The first feature is the general similarity in the observed 
spectra for J/0 radiative decays (left five panels) and 
0(25) radiative decays (right five panels). The yields 
for 0(25) decays are lower. As mentioned earlier, this 
is mainly the consequence of the well-known “13% rule” 
of pQCD which stipulates that the hadronic and radia¬ 
tive decays of 0(25) have nearly a factor seven smaller 


widths compared to those of J/0, a ratio which follows 
that for the leptonic decays. The second feature is the 
expected similarity between the corresponding charged 
and neutral final states, tt+tt" and and K'^K~ 

and KgKg. We find that resonance masses determined 
from charged and neutral spectra are always in agree¬ 
ment within their statistical errors. However, the mass 
resolution and statistics is better for the charged final 
states, as expected. This leads to our decision to make 
resonance mass determinations from the charged decay 
spectra, and to fix them for the final analysis of the neu¬ 
tral decay spectra. The limited event statistics also does 
not allow us to make partial wave analysis of the data. 
One additional observation is that the yields for radia¬ 
tive decays to the 7777 final states, though suggestive of 
the population of 02(1525) and /o(1710), are too small 
to lead to any meaningful results for the branching frac¬ 
tions, and we do not attempt to analyze the 7777 spectra. 

We analyze the spectra of invariant mass of pseu¬ 
doscalar pairs for the radiative decays of J/0 and 0(25) 
in terms of non-interfering relativistic Breit-Wigner res¬ 
onances. The resonance masses, M, the counts, N, and 
the background are kept free. It was not found possible 
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FIG. 9. Two particle invariant mass distributions for J/V', V’(2>S') —>■ 7 PS, PS (PS = K^, Kg, rf) decays, for M(PS, PS) < 

3.0 GeV. 


to keep resonance widths also free, and they were fixed 
at their PDG values jl^. The instrumental resolution 
widths, as given in Table HIl were convolved with the 
PDG widths to obtain the peak shapes used in the fits. 


shape 


BW(M) „ ^ 


M2]2 + [Mr(M)]2 


( 1 ) 


We use the following relativistic Breit-Wigner line- 


where Mq and Pq are the non-relativistic values for the 
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FIG. 10. Illustrating fits to the invariant mass distributions for J ^ 77 r 7 r and yKK. 


mass and width of the resonance, respectively. The mass- 
dependent resonance width T(M) is given by 


r(M) = ro 



f p 

■ B{p) ' 

Uoj 

B{po)_ 


( 2 ) 


where Tq is the non-relativistic width, p and pq are 
the momenta of the resonance’s decay products in its 
center-of-mass frame for masses M and Mq, respec¬ 
tively. B{p) is the L — dependent Blatt-Weisskopf bar¬ 
rier factor, defined as B{p) = 1, l/-\/l -|- {RpY, and 
1 / 7/1 -I- {RpY/?) + {RpY/Q for resonances with L = 0,1, 
and 2, respectively, and R = 3 GeV“^. 

The product branching fractions Bi{J/ip,^{2S) —>■ 
yi?) X S 2 (i? —f PS,PS) are obtained as 


-fR)xB2{R ^ PS,PS) = 


Nr/cr 


N{J/^,^{2S)) 

( 3 ) 

where Nr is the number of events in the resonance R, 
for which the efficiency is tR as given in Table im and 
N{J/ip,ijj{2S)) are the number of J/^/ and '0(25') in the 
data. 


A. Resonances in J/^p Radiative Decays to 
Pseudoscalar Pairs 

The invariant mass spectra for J/ip radiative decays 
and fits to the observed resonances are shown in Fig. IIOI 


Table nni lists the fit results for the resonances. We note 
that in general, the results shown in Table IIIII show that 
for both TTTT and KK the masses determined for charged 
and neutral decays are consistent within their respective 
statistical errors, and the branching fractions are consis¬ 
tent with the isospin expectation of factor two between 
charged and neutral decays. 

We now discuss the results for individual resonances 
separately. 

The scalar fo(500) in J/'0 decay 

The TTTT spectrum: In the tt+tt” spectrum there is a non¬ 
trivial yield in a broad bump in the low energy tail region 
of the strong peak of p(770). It is tempting to ascribe it 
to /o(500). However, such an assignment is unsupport- 
able because we note that there is no evidence for a cor¬ 
responding enhancement in the 7 r° 7 r° spectrum which is 
cleaner in that mass region because of the absence of the 
/ 9 ° peak. 

The vector p{770) in J/ip decay 

As mentioned in Sec. IV the prominent peak for p{77Q) 
arises in the tt+tt” spectrum because of the decay with 
large branching fraction B{J/ip —>■ p° 7 r°) = 5.6(7) x 10“^, 
—>■ TT+TT”, with one photon missing from the tt^ —>■ 
77 decay. As expected, the peak is missing in the tt^tt^ 
spectrum. 

The tensor resonance /2(1270) in J/ip decay 

The TTTT Spectra: This is the strongest excited peak 
in both the tt+tt” and 7 r° 7 r° spectra. Fitted sepa- 
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TABLE III. Fit results for J/^ decays to 7R, R PS, PS. For each observed resonance R, results are presented for decays 
to 7r"'"7r“ {x^/d.o.f = 1 . 36 ), decays to {x^/ d.o.f. = 1 . 00 ), decays to K~ (x^/d.o.f. = 1 . 38 ), and decays to K%K% 

{x^/d.o.f. = 1 . 00 ). For fits to and K~, the resonance masses were kept free. For all fits their widths were fixed at 

their PDG 2014 values. The summed branching fractions are Bi x — BiX B 2 (jv'^n~)+Bi x B2(jv°t^°), and Bi x B 2 {KK) = 

[Bi X B2{K+K-)+Bi X B2{K°sK°s)] x ( 4 / 3 ). 


J/'0 —>■ TR- 
R, 

M(R), r(R), MeV 



R — >■ TTTT 



M(MeV) 

N 

BixBzxlO® 

/2(1270), 

1275(1), 185.0(1:1) 

7r’*"7r 

TTTT 

1259(4) 

1722(61) 

688(36) 

2410(71) 

108.8(39) 

65.6(34) 

174.4(52) 

/o(1500), 

1505(6), 109(7) 

7r'^7r~ 

^00 

TT TT 

TTTT 

1447(16) 

163(36) 

11(18) 

174(40) 

11.0(24) 

1.1(17) 

12.1(29) 

/o(1710), 

1720(6), 135(8) 

7r'''7r~ 

^0 0 

TT TT 

TTTT 

1744(7) 

381(34) 

102(18) 

483(38) 

27.9(25) 

9.3(16) 

37.2(30) 

/o(2100), 

2103(8), 209(19) 

7r’’"7r“ 

TTTT 

2090(10) 

529(39) 

215(23) 

744(45) 

44.3(33) 

18.1(19) 

62.4(48) 

J/'0 —>■ tR- 
R, 

M(R), r(R), MeV 



R ^ KK 



M(MeV) 

N 

Z3ix 132X10® 

/o(1370), 

1350(150), 346(77) 

K+K- 

K%K% 

KK 

1360(31) 

430(90) 

48(15) 

478(91) 

23.1(48) 

8.3(26) 

41.9(73) 

f 2 (1525), 

1525(5), 73(i) 

K+R- 

K%Kl 

KK 

1532(3) 

644(47) 

106(15) 

750(49) 

34.8(25) 

18.4(24) 

70.9(46) 

/o(1710), 

1720(6), 135(8) 

K+K- 

K%K% 

KK 

1706(4) 

1004(47) 

185(18) 

1189(50) 

56.2(26) 

32.0(31) 

117.6(54) 

/o(2200), 

2189(13), 238(50) 

K+K- 

K%K% 

KK 

2206(12) 

381(35) 

109(19) 

490(40) 

26.1(19) 

17.8(31) 

58.6(49) 


rately we obtain for tt+tt decay the mass ) = 

1259(4) MeV, and Bi x B 2 {Tr+TT-) x 10^ = 108.8(39) and 
for 7r°7r° decay and Bi x S2(7r°7r°) x 10® = 65.6(34). The 
total branching fraction for tttt decay is therefore Bi x 
S2(7r7r) x 10® = 174.4(52). These results compare with 
M = 1262(8) MeV, Si xi32(7r+7r-) x 10® = 91.4(7)(148)), 
Bi X S2(7r°7r°) x 10® = 40.0(9)(58), and Si x B 2 {mr) x 
10® = 131.4(11)(159), obtained by BES II [T^ for directly 
produced J/4’- We note that the BES spectra have much 
larger backgrounds which may account for the substan¬ 
tial difference from our results. 

The KK Spectra: The /2(1270) is not strongly excited 
in KK decays. In the K'^K~ spectrum there is a promi¬ 
nent bump in the low energy tail of the strong peak due 
to the excitation of /2(1525) resonance. It can not be fit¬ 
ted with a single peak corresponding to /2(1270), and re¬ 
quires including another larger mass resonance, for which 


the best candidate is the controversial / o ( 1370 ), which 
we discuss below. Because of the isospin difference of 
factor two, and the nearly factor two smaller efficiency, 
the corresponding enhancement the KgKg spectrum has 
too small statistics to make independent analysis. We do 
not fit for /2(1270) —>■ KK, but include its small contri¬ 
bution by appropriately scaling down iV^^(/2(1270)) by 
the PDG O ratio S(/2(1270) ^ KK)/B{f2{l270) 
tttt) = 0 . 054 ( 6 ), to obtain (/ 2 ( 1270 )) = 82 and 

(/ 2 ( 1270 )) = 13 . We do not, however, report 
Bi X B2{KK) for 72(1270). BES II reports that there 
is evidence for / o ( 1270 ) decays into KK, but also does 
not report any numerical results [ 2 ^. 

The scalars /o(1370) and /o(1500) in J/i/> decay 

The TTTT Spectra: The resonance / o ( 1370 ) is considered to 
be the crucial ingredient in the search for scalar glue- 
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balls [4,5]. However, its mere existence is considered 
questionable. As a result, PDG (2014) estimates the 
mass of /o(1370) to be in the range of 1200 — 1500 MeV 
and its width to be between 300 — 500 MeV. In both 
the 7 r''" 7 r“ and spectra a small, relatively narrow, 

“bump” is visible in the high energy tail of the large 
/2(1270) peak. Since the well-known resonance /o(1500) 
is expected in the same tail region, it is difficult to as¬ 
sign the observed bump uniquely to either fn(1370) or 
/o(1500), This problem was also noted by Belle [l^. Our 
attempt to make the best possible choice is based on the 
following exercise. 

We analyze our tt+tt” spectrum assuming that a rel¬ 
atively narrow resonance with width between 100 and 
200 MeV is responsible for the “bump” we observe in 
the high energy tail region of /2(1270). As noted be¬ 
fore, we do not have enough statistics to fit the observed 
structures keeping both the mass and width of this res¬ 
onance free. We therefore assume five different values 
of the width, T = 100,125,150,175, and 200 MeV, and 
fit the 7 r^ 7 r“ spectrum for the best value of the mass 
corresponding to each assumed value of the width. The 
corresponding best fit values are found to be respectively: 
M(MeV) = 1452(13), 1445(15), 1439(17), 1432(19), and 
1424(23). These results indicate that the observed nar¬ 
row “bump” can be only assigned to /o(1500) (for which 
many recent measurements report M r; 1470 MeV [^1. 
and not to /o(1370). This does not mean that a broad 
/o(1370) does not exist, but only that we have no evi¬ 
dence for its existence in our tttt spectra which are dom¬ 
inated by the /2(1270) resonance. 

Our attribution of the “bump” to the excitation of 
/o(1500), assuming the PDG values of the width P = 
109 MeV [Hi, leads to M(/o(1500)) = 1447(16) MeV, 
Bi X X 10^ = 11.0(24) and Bi x B2{tt°tt°) x 

10^ = 1.1(17). The summed branching fraction is 
Bi X B 2 {t^t^) X 10^ = 12.1(29). BES II has analyzed 
their data for radiative decays of directly produced J/^/> 
into pion pairs, and reported M(1500) = 1466(21) MeV, 
Bi X H 2 ( 7 r+ 7 r-) x 10^ = 6.7(2)(30), Bi x B 2 {tt°tt°) x 10^ = 
3.4(3)(I5), and Bi x B 2 {'ktt) x 10^ = I0.1(4)(34) [Hj. 

The KK Spectra: Both the K^K~ and KgKg spec¬ 
tra are dominated by the /2(1525) and /o(I7I0) reso¬ 
nances. However, there is appreciable yield spread over 
more than 200 MeV in the K'^K~ spectrum in the region 
of the low energy tail of /2(1525). To determine whether 
this yield should be assigned to the resonance /2(1270) 
or /o(1370), we have done the same exercise as described 
above for /o(1500). We fit the K^K~ spectrum to deter¬ 
mine values of masses in this region with assumed values 
of width, P = 200, 250, 300, 350, and 400 MeV. The best 
fit masses are found to vary from 1330 to 1353 MeV. 
Allowing the width to be a free parameter results in 
M = 1360(31) MeV, P = 346(77) MeV. We therefore 
assign this enhancement to /o(1370). The K'^K~ spec¬ 
trum leads to Bi x B 2 {K^K~) x 10^ = 23.1(48), and the 
K'gK^ spectrum leads to BixB 2 {K^K^) x 10^ = 8.3(26). 
The summed result is Bi x B 2 {KK) x 10® = 41.9(73), 


taking account of a factor 4/3 for the undetected AT® AT® 
decay, as noted for all kaon decays in the Table [Till cap¬ 
tion. BES H also observes kaon yields in this general mass 
region, but does not quote any numerical results [T^ . 

The tensor /2(1525) in J/i/j decay 

There is no evidence for the excitation of /2(1525) in 
the TTTT decay spectra. In the KK spectra /2(1525) is 
prominently excited, with the K~^K~ spectrum lead¬ 
ing to M(/^(1525)) = 1532(3) MeV. The K+K- spec¬ 
trum leads to Bi x B 2 {K'^K~) x 10® = 34.8(25), and 
the KgKg spectrum leads to Bi x B2{KgKg) x 10 ® = 
18.4(24). The summed result is Bi x 132 (ATAT) x 10® = 
70.9(46). The result from BES H with large background 
in the invariant mass spectra is Bi x B 2 {KK) x 10 ® = 
34.2(15)(i)(i®®) 0. 

The scalar /o(1710) in J/-!/) decay 

The scalar /o(1710) is prominently excited in both tttt 
decays and KK decays iFig. ITOl) . The 7r'''7r“ spectrum 
leads to the mass M(/o(1710)) = 1744(7) MeV and 
Bi X ^ 2 ( 71 + 71 “) X 10® = 27.9(25), and the 7r°7r° spec¬ 
trum leads to Bi x H 2 ( 7 r® 7 r®) x 10® = 9.3(16), with the 
sum Bi X ^ 2 ( 7777 ) X 10® = 37.2(30). BES H has reported 
Bi X -B 2 ( 7 r+ 7 r-) x 10® = 26.4(4)(75) and Bi x B 2 {'k°tt^) x 
10® = 13.3(5)(88), and therefore Bi x ^ 2 ( 7171 ) x 10® = 
39.7(6)(116) [H. 

The KK Spectra: The K^K~ spectrum leads to 
M(1710) = 1706(4) MeV, and Bi x B 2 {K+K-) x 10® = 
56.2(26), the KgKg spectrum leads to Bi x B2{KgKg) x 
10® = 32.0(31). Our summed result is Bi x B 2 {KK) x 
10® = 117.6(54). The corresponding result of BES II is 
Hi X B2{KK) X 10® = 96.2(29)(2ii)( 281) [l^. 

Of the seven resonances listed in Table m unambigu¬ 
ous results for both tttt and KK decays are observed for 
only one resonance, /o(1710). This is largely due to the 
fact that all the resonances have widths much larger than 
the experimental resolution widths, and they strongly 
overlap. As a result, as described in Sec. IX, only for 
/o(1710) we are able to test the flavor independent char¬ 
acter attributed to pure glueballs. 

Both K^K~ and KgKg spectra are very well defined 
in their high energy tail region, and show no evidence 
for enhancement in the M = 1790 region which has been 
occasionally reported. 

The M — 1800 — 2500 MeV region in Jj-ip decay 

The mass region above /o(1710) is complicated with sev¬ 
eral /o, and /2 resonances reported [26j |. In particular, 
radiative decays of J/ijj have been reported to /o(1800), 
/o(2200), /2(1910), /2(1950), and to the much talked 
about tensor / 2 ( 2220 ). In the tttt and KK mass spec¬ 
tra in the present investigation a broad enhancement at 
M r; 2100 MeV is seen in the tttt decays, and another 
broad enhancement at M r: 2200 MeV is seen in the KK 
decays. The background is generally much larger than in 
the lower mass regions. This makes resonance assign¬ 
ments in this mass region rather speculative. With this 
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FIG. 11. Illustrating fits to the invariant mass distributions for '4>{‘2S) —^ 'yiin and 'yKK. 



FIG. 12. Data for 7r’'"7r“ and K^K~ mass distributions with 
superposed continuum contributions (in blue) estimated by 
normalizing continuum yields observed in off-resonance data 
taken at ^/s Ri 3871 MeV. 


caveat we fit the tttt spectra in this region with /o(2100) 
with the PDG width of 209 MeV, and the KK spectra 
with /o(2200) with the PDG width of 238 MeV. Fits to 
the TTTT spectra lead to M(/o(2100)) = 2090(10) MeV, 
Bi X B 2 {Tr+Tr-) x 10^ = 44.3(33), Bi x B2{Tr°7r°) x 10® = 
18.1(19), and the sum Bi x B 2 {t^tt) x 10® = 62.4(48). Fits 
to the KK spectra lead to M(/o(2200)) = 2206(12) MeV, 
Bi X B 2 {K+K-) X 10® = 26.1(19), and Bi x B 2 {K%K%) x 
10® = 17.8(51), and the total Bi x B 2 {KK) x 10® = 
58.6(49). 


B. Resonances in 'i/'(2S) Radiative Decays to 
Pseudoscalar Pairs 

While searches for glueballs in radiative decays of J/^/> 
have been reported before, similar searches in 1^(23) 
radiative decays have not been reported before. The 
only reports for '0(25') radiative decays which exist 
are a BES II publication with a sample of 4.02 mil¬ 
lion '0(2<S') [id, and an unpublished report with 14 mil¬ 
lion '!/'(2<5') [^. We present in this paper our results for 
the radiative decays observed in a sample of 24.5 mil¬ 
lion ijj{2S). 

As pointed out in Sections IV and V, and shown in 
Fig. 7, the results for '0(25') decays are very similar to 
those of J/0 decays. The event statistics is smaller, due 
primarily to the “13% rule”, and the backgrounds are 
relatively larger. The fits to the invariant mass spectra, 
and the determination of the product branching fractions 
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TABLE IV. Fit results for 'i/>(2S) decays to 7 R, R —>■ PS, PS. For each observed resonance R, results are presented for 
decays to 7 r"'" 7 r“ {x^fd.o.f. = 1.24), decays to (x^ fd.o.f. = 1.19), decays to K~ {x^/d.o.f. = 1.29), and decays to 

KgKg {x^ld.o.f. — 0.94. For fits to 7 r''’ 7 r“ and K~, the resonance masses were kept free. For all fits their widths were 

fixed at their PDG2014 values. The summed branching fractions are Bi x B 2 (tvtv) = x B 2 {Tv'^Tr~) + x l? 2 ( 7 r° 7 r°) and 
Bi X B2{KK) = [Bi X B2{K+K-) +BiX B2{K°sK°s)] x (4/3). 


0(25) -)■ 7 R 

R, 

M(R), r(R), MeV 



R — >■ TTTT 



M(MeV) 

N 

5ix 52X10® 

02(1270), 

7r’*"7r 

1267(4) 

1509(63) 

17.0(7) 

1275(1), 185.0(1:1) 



415(30) 

6.9(5) 


TTTT 


1924(70) 

23.9(9) 

00 (1500), 

7r’’'7r~ 

1442(9) 

261(40) 

3.0(5) 

1505(6), 109(7) 

^0^0 

TT TT 


13(17) 

0.2(3) 


TTTT 


274(43) 

3.2(6) 

00(1710), 

7r’’'7r~ 

1705(11) 

237(30) 

2.8(4) 

1720(6), 135(8) 

^0^0 

TT TT 


53(14) 

0 . 8 ( 2 ) 


TTTT 


290(33) 

3.6(4) 

00 ( 2100 ), 

7r’’"7r“ 

2099(17) 

283(35) 

3.4(4) 

2103(8), 209(19) 



90(17) 

1.4(3) 


TTTT 


373(39) 

4.8(5) 

0(25) -)■ 7 R 

R, 

M(R), r(R), MeV 



R ^ KK 



M(MeV) 

N 

5ix 52X10® 

0o(137O), 

K+K- 

1350(48) 

168(50) 

1.9(6) 

1350(150), 346(77) 

K%K% 


7(12) 

0.4(4) 


KK 


175(51) 

3.1(10) 

02(1525), 

K+K- 

1552(9) 

113(24) 

1.3(3) 

1525(5), 73(i) 

K%Kl 


23(8) 

0.9(3) 


KK 


136(25) 

2.9(6) 

0o(171O), 

K+K- 

1690(8) 

349(30) 

4.1(3) 

1720(6), 135(8) 

K%K% 


26(8) 

0.9(3) 


KK 


375(31) 

6.7(6) 

00 ( 2200 ), 

K+K- 

2188(17) 

203(30) 

2.3(3) 

2189(11), 238(50) 

K%K% 


4(9) 

0.1(3) 


KK 


207(31) 

3.2(6) 


for the ip{‘^S) decays was done in the same manner as for 
J/'0 decays. The results are tabulated in Table ITVl and 
the fits to the charged and neutral decays for tttt and KK 
spectra are shown in Fig. [TT] 

The vectors p(770) and 0(1020) in xp{2S) decay 

In the TT+TT” and K~^K~ invariant mass spectra we note 
the presence of the p(770) and 0(1020) peaks as also indi¬ 
cated by the diagonal enhancements in the corresponding 
Dalitz plots of Fig. [8] These arise due to the production 
of p(770) and 0(1020) vector mesons via initial state ra¬ 
diation. To confirm this, we have analyzed our off-'0(2iS') 
resonance data at y/s = 3671 MeV with the integrated 
luminosity of 20.6 pb“^. To compare the 7r''"7r“ and 
K'^K~ yields observed in these data with those observed 
in 0(25) data we multiply the yields by the relative lumi¬ 
nosity factor £(0(25))/£(3671) = 53/20.6 = 2.57, and 


the factor (s'/s)^ = (3686/3671)^ = 1.06. The normal¬ 
ized yields are shown in Fig. [12] superposed on the 0(25) 
data. It is seen that the normalized contributions fully 
account for the yield of /o(770) in the tt+tt” spectrum, 
and for the yield of 0(1020) in the K'^K~ spectrum. 
It is also seen that the background contribution due to 
this ISR yield is predicted to be negligible for masses 
> 1 GeV. 

The tensor /2(1270) in 0(25) decay 

As in J/0 decay, the 02(1270) is prominently excited 
in both 7r+7r“ and 7r°7r° decays. On the high energy 
shoulders in both spectra we again see enhancement 
which is well fitted with the narrow /o(1500) rather 
than narrow or wide /o(1370). The fit to 02(1270) re¬ 
sults in M = 1267(4) MeV, Bi x 52(7r+7r-) x lO'^ = 
17.0(7), Bi X 52 (7r°7r°) x 10^ = 6.9(5), and the sum 
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Bi X SzCtttt) X lO'^ = 23.9(9). 

In the K~^K~ spectrum there is a broad enhancement 
below 1500 MeV which cannot be assigned to /2(1270). 
Instead, it is well fitted with a wide resonance with mass 
around 1350 MeV, as discussed below. 

The scalar /o(1370) in 'ip{2S) decay 

In the J/V' decay to tttt it was not possible to discern the 
possible excitation of /o(1370) because of the strong ex¬ 
citation of /2(1270). As described in Sec. V.A, the small 
shoulder seen on the high energy tail of /2(1270) could 
be assigned to /o(1500) but not to /o(1370). In the KK 
spectra for 4’i‘^S) the /2(1270) excitation is much weaker. 
This allows us the opportunity to fit the mass region be¬ 
low 1525 MeV. The 'tjj(2S) spectra have much smaller 
statistics than the J/'i/f spectra, which makes it difficult 
to determine the properties of a wide, weakly excited res¬ 
onance like the /o(1370). To explore its possible excita¬ 
tion, we have fitted the K'^K~ spectrum with the width 
of /o(1370) fixed to T = 250,300,350 MeV. The result¬ 
ing masses were M = 1356(41), 1353(45), 1349(48) MeV, 
i.e., in essential agreement with the determination of 
M = 1360(31) MeV and F = 346 MeV for J/ip decay. 

For our final fit, we therefore assume a state with width 
of 346 MeV. The fit to the K'^K~ spectrum yields a mass 
oi M = 1350(48) MeV. We identify it with /o(1370). 
The fit results to Bi x B 2 {K'^K~) x 10® = 1.9(6), Bi x 
B2{K%K%) X 10® = 0.4(4), and the total Bi x B2{KK) x 
10® = 3.1(10). This fit has been made excluding the two 
bins near M ^ 1350 MeV, which would correspond to a 
very narrow state of F <50 MeV, for which no plausible 
candidate exists. 

The scalar /o(1500) in if){2S) decay 

The absence of /2(1525) in tttt decay allows us to fit for 
/o(1500) with the results: M(/o(1500)) = 1442(9) MeV, 
Bi X ,B2(7r+7r-) x 10® = 3.0(5), Bi x B2(7r°7r°) x 10® = 
0.2(3), and the sum Bi x B 2 {tt'k) x 10® = 3.2(6). 

The strong excitation of /2(1525) in KK decays makes 
it impossible to say anything about possible excitation of 
/o(1500) in KK decays. 

The tensor /2(1525) in 'ip(2S) decay 

No evidence for the tttt decay of /2(1525) was found in the 
radiative decay of and none is found in the radia¬ 
tive decay tp{2S) —>■ 77 r 7 r, either. In K^K~ decay strong 
excitation of /2(1525) is observed, and the fit to the 
K'^K~ spectrum leads to M(/2(1525)) = 1552(9) MeV 
andSiX,B 2 (A:+A:-)xlO® = 1.3(3). The fit to the 
spectrum leads to Bi x B 2 {KgKg) x 10® = 0.9(3). The 
total branching fraction is Bi x B 2 {KK) x 10® = 2.9(6). 

The scalar /o(1710) in xp{2S) decay 

The scalar /o(1710) is strongly excited in both tttt 
and KK decays. Fits to the tttt spectra lead to 
M(/o(1710)) = 1705(11) MeV, Bi x ,B2(7r+7r-) x 10® = 
2.8(4), and Bi x i32(7r°7r®) x 10® = 0.8(2), and the sum 
Bi X B 2 (t^t^) X 10® = 3.6(4). Fits to the KK spectra lead 


to M(/o(1710)) = 1690(8) MeV, BixB 2 {K+K-) x 10® = 
4.1(3), Bi X B 2 {K%K%) X 10® = 0.9(3), and the total 
Bi X B2{KK) X 10® = 6.7(6). 

The M — 1800 — 2500 MeV region in xp{2S) decay 

As mentioned for J/ij} decays, in this region backgrounds 
are large and only a large, wide peak is observed in both 
TTTT and KK spectra. Each of them may be composite 
of the excitation of more than one resonance. If ana¬ 
lyzed as single resonances, fits to the tttt spectra lead to 
M(/o(2100)) = 2099(17) MeV, Bi x B2 {'k+tt-) x 10® = 
3.4(4), Bi X ,B2(7’'°7r°) x 10® = 1.4(3), and the sum 
Bi X B 2 {Tnr) X 10® = 4.8(5). Anisovich and Uman have 
reported the observation of a resonance with this mass 
in pp annihilations [1^. Fits to the KK spectra lead to 
A/(/o(2200)) = 2188(17) MeV, Bi x B2{K+K-) x 10® = 
2.3(3), Bi X B 2 {K%KI) X 10® = 0.1(3), and the total 
Bi X B 2 (KK) X 10® = 3.2(6). Augustin [l^ has reported 
excitation of a resonance with this mass in radiative de¬ 
cay of J/'0. 


C. The ^(2230) Tensor Glueball Candidate 


While the main interest in glueball searches has been 
in the scalars, special interest has existed in the claim 
for a narrow tensor state with mass ^ 2230 MeV, and 
width F Ri 20 — 30 MeV. This claim has a check¬ 
ered history. Originally claimed by Mark III [T^ in 
the radiative decay of J/4’, and the decay of ^(2230) 
into K'^K~ and KgK^, it has gone through several 
claimed observations fl5 . IT^ . . and a number of non¬ 
observations [i3, mi, [iMlf . The latest of these consist 
of its observation by ITEP in 7r~p formation and decay 
into KgKg [ 2 ^, and its non-observation by BaBar in the 
radiative decay of J/ijj into AT+AT” and KgKg [^, and 
its non-observation by BES III in the radiative decay 
of J/tjj and decay into rjrj mi- The present measure¬ 
ment allow us to search for it in the radiative decays 
of both J/'f/j and tf}{2S), with its decays into tttt, KK, 
and 7777 pseudoscalar pairs. In the analyses of the invari¬ 
ant mass spectra as described above, we add a resonance 
with M = 2230 MeV, and widths of 20 and 50 MeV, and 
determine upper limits at 90% confidence level for the 
product branching fractions for the existence of ^(2230). 
The results are listed in Table |Vl They can be summa¬ 
rized as upper limits of < 5 x 10“® in J/ip decays and 
< 6 X 10“® in 7 / 1 ( 25 ') decays. Our limits in J/t/j decays 
are consistent with < 2.3 x 10“® reported by DM2 
for K~^K~ and KgKg decays, and < 4 x 10“® reported 
by BaBar m| in K^K~ and KgKg decays. No evidence 
for its existence was found by BES III in radiative decay 
of a very large sample of 2.25 x 10® J/t/; to 7777 mi- 
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TABLE V. Upper limits at 90% confidence level for the production of the presumed tensor glueball candidate ^(2230) in the 
reactions J/'ip,'tp{2S) —>■ 7 R, R —>■ PS, PS. 


r(C), MeV 

90% CL UL 

N B 1 XB 2 

xlO® 

Bi X B 2 (ref.) 
xlO® 

J/ijj —>■ 7(7r+7r~) 20/50 

-)■ 7(7r°7r°) 20/50 

7 (A'+A-) 20/50 

^ 7(^s^s) 20/50 

1/(25') —>■ 7(7r^7r~) 20/50 

7(7r°7r'’) 20/50 

-)■ 'y(K+K-) 20/50 

^ 'y{K°K°) 20/50 

23.0/45.1 2.6/5.2 
16.0/23.2 1.3/1.9 

25.4/44.6 1.7/3.1 

11.5/19.1 1.2/2.0 

22.9/31.1 0.32/0.43 
15.7/23.6 0.26/0.40 
18.9/38.0 0.21/0.43 
10.7/15.6 0.37/0.55 

5.6(^^) [ 1 ^ 

3.3(11) [18J, 4.2(«) [ 1 ^ 
2.3 (95% CL) [17] 
3.1(15) [18], 2.7(1!) [1^ 
1.6 (95% CL) [17] 


TABLE VI. Summary of results for J/'i/) and ip{2S) radiative decays to light quark resonances. For entries with two uncertainties, 
the first listed uncertainties are statistical and the second are systematic as described in Sec. VIII and listed in Tables VIII 
and IX. 



From J/ip 

M (MeV) B 1 XB 2 X 10® 

From 1/(25) 

M (MeV) B 1 XB 2 X 10® 

62(1/(25)) 

B 2 (J/i/) 

(%) 

/2(1270) ->■ TTTT 

1259(4)(4) 

174.4(52) (122) 

1267(4) (3) 

23.9(9)(9) 

13.7(7)(11) 

/o(1370) ->■ KK 

1360(31)(28) 

41.9(73)(134) 

1350(48)(15) 

3.1(10)(14) 

7.4(27)(41) 

/o(1500) — >■ TTTT 

1447(16) (13) 

12.1(29)(24) 

1442(9)(4) 

3.2(6)(2) 

26.4(80)(55) 

/!(1525) ->■ KK 

1532(3)(6) 

70.9(46)(67) 

1557(9) (3) 

2.9(6)(3) 

4.1(9)(6) 

/o(1710) ->■ TTTT 

1744(7)(5) 

37.2(30)(43) 

1705(11)(5) 

3.6(4)(5) 

9.7(13)(18) 

KK 

1706(4)(5) 

117.6(54)(94) 

1690(8)(3) 

6.7(6)(6) 

5.7(6)(7) 

/o(2100) ->■ TTTT 

2090(10)(6) 

62.4(48)(87) 

2099(17) (8) 

4.8(5)(9) 

7.7(10)(18) 

/o(2200) KK 

2206(12)(8) 

58.6(49)(120) 

2188(17)(16) 

3.2(6)(8) 

5.5(11)(18) 


TABLE VII. Product branching fractions for Xc 0,2 —>■ PS,PS 
decays (based on Ref. [s^). 


B 1 XB 2 X 10® 

XcO 

Xc 2 

XcJ 

—>■ TT^TT 

58.8(30) 

14.8(8) 

XcJ 

^ TT TT 

27.1(30) 

6.3(7) 

XcJ 

TTTT 

85.9(42) 

21.1(11) 

XcJ 

->■ K+K- 

58.5(32) 

10.6(6) 

XcJ 

^K%K% 

32.2(18) 

4.9(4) 

XcJ 

KK 

120.9(49) 

20.7(10) 

XcJ 

VV 

29.2(31) 

4.8(7) 

XcJ 

v'v' 

19.4(20) 

0.53(31) 

XcJ 

m' 

1.46(56) 

0.13(29) 


VI. SUMMARY AND COMPARISON OF 
RESULTS FOR RADIATIVE DECAYS OF J/V> 
AND il){2S) TO RESONANCES WITH MASSES 

M < 2500 MEV 

In Table EH we compare the results for all resonances 
observed in the radiative decays of J/r/i and ^/’(2S') to 
pairs of charged and neutral pseudoscalars. Several fea¬ 
tures are notable in Table ED 

As mentioned in Sec. V.B, the results for J/ip and 
il’{2S) decays are very similar for both charged and neu¬ 
tral decays. The masses are generally in agreement well 
within their statistical errors. Branching fraction results 
from BES, when available, are generally not in good 
agreement with ours. The BES results are based on di¬ 
rectly produced J/V', and as mentioned before, the BES 
invariant mass spectra contain large backgrounds due to 
ISR contributions. 

The So-called 13% Rule: The “13% rule” mentioned 
earlier is generally invoked for hadronic decays of vec- 
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FIG. 13. Two particle invariant mass distributions for se¬ 
lected V'(2S') ^ 7 PS,PS, (PS = 7v^,Tr°,K^,K°s,v) decays, 
for the region M > 2.5 GeV. The XcO and Xc 2 states are 
clearly seen. No evidence for any other states is seen. 


tor resonances of charmonium J/tp, ip(2S), ... ip(nS) be¬ 
cause both their leptonic decays via one photon and their 
hadronic decays via three gluons are proportional to the 
wave function at origin. Actually, the radiative decays 
to hadrons, via a photon plus two gluons have the same 
proportionality to the wave function at origin, so that 
we expect, for example, the ratio of radiative decays to 
pseudoscalar pairs 

B{ip{2S) 7PS,PS)/S(J/V’ ^ 7PS,PS) 

= B(V'(2S') ^ e+e-)IB{Jlip e+e”) 

= 0.789 X 10-V5-971 x IQ-^ = 13.44(30)% 

In Table EJ we list this ratio for the various resonances 
we observe. The ratio is found to vary between 3.6(8)% 
and 26.2(84)%. This is similar to the variation which has 
been observed for the pure hadronic decays [s^- 


VII. RESONANCES WITH MASSES > 2500 MEV 


A. Search for Glueballs with Masses in the 
2.5 — 3.3 GeV Region 

One of the main reasons for our studying ip{2S) radia¬ 
tive decays in addition to those oi J/ip was to extend our 
glueball searches to masses in the region 2.5 — 3.3 GeV. 
As stated earlier, and illustrated in Fig.[Tl lattice calcula¬ 
tions predict a = 0 glueball with mass M{Q '') = 
2560(35)(120) MeV, and a = 2 '' glueball with 
mass M(2-+) = 3040(40)(150) MeV. These 0"+ and 

2 ^ positive charge conjugation states can be excited 
in radiative decays of 'tp{2S), and we have searched for 
them and any other enhancements in the entire region, 
M = 2.5 — 3.3 GeV, in their decays into pseudoscalar 
pairs, assuming three different values of their widths, 
r(A:) = 50 MeV, 100 MeV and 150 MeV. The results are 
displayed in Fig. [TH No notable enhancement is found 
anywhere in the mass region, and upper limits for their 
branching fractions at 90% confidence level have been 
determined. The upper limits for all product branching 
fractions Bi x B 2 are determined to be all < 1 x 10 “®, 
except for decay to 7777 for which they become as large as 

3 X 10“® because of the very small number of observed 
counts and small B{r] —>■ 77 ) ~ 0.4. 


B. The XcO and Xc 2 resonances of charmonium 


As shown in the Dalitz plot of Fig. [H and invari¬ 
ant mass projections in Fig. [T3] for 'ip{2S) radiative de¬ 
cay, we observe strongly populated XcO and Xc 2 states 
of charmonium in their decay info tttt, KK, and rjrj. 
The branching fractions for these decays were reported 
by us earlier in the GLEO publication of Ref. [ 3 ^ . 
The S 2 (XcO,c 2 —>■ PS,PS) branching fractions reported 
there, multiplied by Bi{'ip{2S) —>■ 7Xco) = 9.22% and 





























17 



' 2.5 2.6 2.7 2.8 2.9 




FIG. 14. Upper limits at 90% confidence level of R) x B 2 {R ^ PS,PS) as fnnction of resonance mass, M. The 

differently colored cnrves, from top to bottom, correspond to upper limits for assumed widths r(R) = 150, 100 and 50 MeV. 
The two-pseudoscalar invariant mass spectra are shown as the shaded histograms, with the scale of the observed connts on the 
right ordinate. The thick horizontal lines illnstrate the mass ranges predicated in Ref. [13] for the glueball candidates with 

jPC ^ Q-+ 


B2{ip{‘2.S) —>■ 7 Xc 2 ) = 9.33% used there, are quoted as 
B 1 XB 2 X 10^ in Table VII. 

Ochs @ has emphasized that while it is difficult to 
study the flavor-independent decay of two-gluon glueball 
candidates to pseudoscalar pairs, the corresponding de¬ 


cays of Xco{J = 0 ) and Xc 2 {J = 2 ), which are mediated 
by two gluons, provide an excellent opportunity to test 
the flavor independence ansatz. The values for tttt, KK, 
and 7777 decays listed in Table VII lead to the ratios below: 
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TABLE VIII. Summary of systematic uncertainties on branching fractions (in units of percent). 



J/tp -fTV 

/yTT^TT^ 


tp{2S) —>■ 77r+ 

1 

C 

TT*' 


/2(1270) /o(1500) /o(inO) /o(2100) 

/2(1270) /o(1500) /o(1710) /o(2100) 

Reconstruction & PID 

5/5 

5/5 

5/5 

5/5 

4/5 

4/5 

4/5 

4/5 

Kinematic Fit 

1/2 

10/6 

1/8 

A/12 

1/1 

5/29 

3/3 

23/2 

B-W Widths 

1/1 

6/8 

3/3 

7/5 

1/1 

4/10 

10/4 

7/8 

Background 

1/1 

3/13 

6/5 

7/6 

1/2 

2/2 

5/9 

10/11 

Binning 

1/1 

2/5 

2/1 

2/2 

1/1 

1/29 

2/1 

2/1 

pn Rejection 

5/- 

16/- 

8 /- 

5/- 

— 

— 

— 

— 

Number of ipi^S) 

2/2 

2/2 

2/2 

2/2 

2/2 

2/2 

2/1 

2/2 

Total 

8/6 

21/18 

12/11 

13/15 

5/6 

8/43 

13/12 

26/15 


JIP, 


K-hK'^sKl 

^P{2S)^xK^K-/ xKlKl 


/o(1370) /2(1525) /o(1710) /o(2200) 

/o(1370) /2(1525) /o(1710) /o(2200) 

Reconstruction & PID 

6/5 

6/5 

6/5 

6/5 

5/5 

5/5 

5/5 

5/5 

Kinematic Fit 

28/3 

4/5 

2/3 

5/2 

10/7 

8/4 

2/10 

1/30 

B-W Widths 

20/11 

5/4 

3/3 

8/16 

35/24 

7/6 

4/4 

13/3 

Background 

11/6 

1/1 

1/2 

2/11 

26/100 

6/5 

7/10 

22/100 

Binning 

3/4 

1/3 

1/1 

1/2 

24/15 

6/1 

3/2 

2/40 

K * K Rejection 

28/- 

5/- 

6 /- 

18/- 

— 

— 

— 

— 

Number of ip{2S) 

2/2 

2/2 

2/2 

2/2 

2/2 

2/2 

2/2 

2/2 

Total 

46/18 

10/9 

9/7 

21/20 

51/104 

15/10 

10/16 

26/112 


Pseudoscalar pairs 
Two-gluon decays (theoretical) 

Xco decays 
Xc 2 decays 

XcO decays (divided by phase space factor p) 
Xc 2 decays (divided by phase space factor p^) 


The branching fraction ratios for XcO decays are seen to 
agree within experimental uncertainties with the theoret¬ 
ical expectations for 2-gluon decays. The Xc 2 ratios also 
agree with these expectations, except for decays to rj'r]'. 


VIII. SYSTEMATIC UNCERTAINTIES 

Our estimates of systematic uncertainties in product 
branching fractions J/ij) ^ yTTTr, ^KK, and iIj{2S) —>■ 
•yiriT, jKK, are listed in Table Ivnil We estimate uncer¬ 
tainties due to reconstruction and particle identification 
as follows. We find a reconstruction uncertainty of 1% 
per charged particle, with an additional 1% and 2% per 
identified and K^, respectively, for the use of dE/dx 
and RICH information. We find a —)> tt+tt” iden¬ 

tification uncertainty of 1%. We also find uncertainties 
of 2% per 7r° and 2% for the radiative photon. We add 
the contributions for each mode in quadrature to obtain 
the systematic uncertainties for this source. We estimate 


TTTT 

KK 

m 

v'v' 

m' 

= 1 

1.33 

0.33 

0.33 

0 

= 1.00(5) 

1.41(6) 

0.34(4) 

0.23(3) 

0.02(1) 

= 1.00(5) 

0.98(6) 

0.23(4) 

0.025(15) 

0.006(15) 

= 1.00(5) 

1.47(9) 

0.36(4) 

0.27(3) 

0.02(1) 

= 1.00(5) 

1.18(8) 

0.28(4) 

0.058(34) 

0.010(22) 


the uncertainty in our requirement that the kinematic fit 
have x^ < 25 by analyzing data with the requirements 
X^ < 20 and x^ < 30 j and we take the largest varia¬ 
tions as the systematic uncertainty due to this source. 
In resonance fits, our fixed parameters are the Breit- 
Wigner widths of the resonances and the MC-determined 
mass resolution. To estimate the uncertainties due to 
the fixed resonance widths, we refit the data, varying the 
widths by ±lcr of their PDG values [2^, and we take 
the largest variation for each resonance as the system¬ 
atic uncertainty due to this source. We also refit the 
data with the MC-determined mass resolution varied by 
±10%, and find that this has no measurable effect our re¬ 
sults, since the resonances in our fits are all much wider 
than the mass resolutions. The uncertainty due to the 
backgrounds in the fits is obtained by varying background 
normalization by ±lcr around the central value obtained 
from the fits. The uncertainty due to the 25 MeV binning 
of the invariant mass spectra is obtained by analyzing the 
data with 15 MeV and 5 MeV bins and taking the largest 
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TABLE IX. Summary of systematic uncertainties on resonance masses (in units of MeV). 




J/0- 

->■ 777^77 



0(25) - 

->■ 777''' 77 



/2(1270) /o(1500) /o(1710) /o(2100) 

/2(1270) /o(1500) /o(1710) /o(2100) 

Mass Calibration 

2 

2 

2 

2 

2 

2 

2 

2 

Kinematic Fit 

1 

7 

2 

5 

1 

1 

3 

5 

B-W Widths 

1 

2 

1 

1 

1 

2 

2 

3 

Backgronnd 

1 

2 

1 

1 

1 

1 

1 

4 

pn Rejection 

3 

9 

4 

1 

— 

— 

— 

— 

Binning 

1 

5 

1 

3 

1 

2 

2 

3 

Total 

4 

13 

5 

6 

3 

4 

5 

8 



J/0- 

> -yK+K- 



0(25) ->■ 7A+A- 



/o(1370) /2(1525) /o( 1710) /o(2200) 

/o(1370) /2(1525) /o( 1710) /o(2200) 

Mass Calibration 

2 

2 

2 

2 

2 

2 

2 

2 

Kinematic Fit 

2 

1 

1 

3 

10 

2 

2 

3 

B-W Widths 

25 

2 

1 

3 

8 

1 

1 

3 

Backgronnd 

5 

1 

1 

1 

5 

1 

1 

5 

K*K Rejection 

9 

4 

4 

6 

— 

— 

— 

— 

Binning 

7 

2 

1 

1 

5 

1 

1 

3 

Total 

28 

6 

5 

8 

24 

3 

3 

6 


variation as the systematic uncertainty due to this source. 
The uncertainty in the number of 'tp(2S) is 2%. 

We estimate the uncertainty in our mass calibration by 
fitting the XcO and Xc 2 resonances observed in the 'ip{2S) 
data. We find that the XcO and Xr.i masses in all decays 
agree with their known values 2^ within ±2 MeV, and 
we take ±2 MeV as the systematic uncertainty in our 
mass calibration. Other contributions to uncertainties 
in mass determinations from the sources described above 
are also listed in Table El 

The sums in quadrature of all systematic uncertainties 
are also listed in Tables IVIIII and El For the upper limits 
on product branching fractions reported in Table El we 
take the largest upper limits obtained from these varia¬ 
tions as our hnal upper limit results. 


invariant mass spectra for J/tj) radiative decay also con¬ 
tributes to the differences, particularly for weakly excited 
states. Our spectra for J/tj; radiative decay are free from 
ISR background contributions. 

A resonance with mass ~ 1800 MeV has been sug¬ 
gested in hadronic decays of J/\j} 01 • We do not find 
evidence for such a state in either J/ijj oi ip{2S) radia¬ 
tive decays. 

Our results for '4’{2S) and J/ij) decays are generally 
in agreement. The ratios of '0(25')/J/'0 branching frac¬ 
tions vary between ^ 5% to 15%. Similar variations from 
the 13% pQCD rule have been observed before for other 
decays [^ . 


B. Results Relating to the Glueballs 


IX. DISCUSSION OF RESULTS AND 
CONCLUSIONS 

A. Branching Fractions for Light Quark 
Resonances 

Our final results for light quark scalar and tensor reso¬ 
nances in the mass region 1100 — 2200 MeV are summa¬ 
rized in Table EH As mentioned earlier, for J/0 radiative 
decays, results from BES measurements often differ sig- 
nihcantly from our results. The differences in the results 
can be attributed mainly to the intrinsic differences be¬ 
tween the methods of analysis. Our analyses are based 
on Breit-Wigner resonance formalism. BES results are 
based on partial-wave analyses, which leads to large un¬ 
certainties depending on the number of resonances used 
in the analysis. The large ISR backgrounds in the BES 


1. The 2++ Glueball ^(2230) 

We find no evidence for the claimed narrow tensor state 
^(2230), and establish limits for its existence for assumed 
widths of 20 and 50 MeV in the different decay modes 
in terms of the product branching fraction at 90% confi¬ 
dence level of -B( J/0 —>■ 7 ^ —>■ 7 (PS,PS)) < 5 x 10“® and 
.6(0(25) ^ 7 ^ ^ 7 (PS,PS)) < 0.6 x 10-0 

2. The 0 ^ and 1 ^ Glueballs 

As shown in Fig. 1141 we do not find any resonance 
enhancement in tttt, KK, or 7777 decay in radiative decays 
0(25) in the mass region 2.5 —3.3 GeV in which 0 ^ and 
1 ^ glueballs are predicted by lattice calculations [l^ . 
At 90% confidence level, upper limits for the product 
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branching fractions between 10 ® and 10 ® are obtained 
for assumed widths 50 — 150 MeV. 


3. The Scalar /o(1370) 

We find convincing evidence for the existence of 
/o(1370), although we are able to identify it only in its 
decay to KK. Our inability to observe it in its tttt decay 
is due to the strong excitation of /2(1270) in its vicinity. 
The existence of /o(1370) has been considered to be the 
crucial evidence for the supernumery nature of scalars in 
this mass region, and therefore of the existence of a scalar 
glueball, albeit mixed with other scalars. 


4- Flavor Blind Decays 

As mentioned in the introduction, the relative decay 
rates to pseudoscalar pairs tttt : KK : rjrj : rj'r]' of 3 : 
4 : 1 : 1 is considered to be the defining signature of a 
pure, unmixed glueball. Our measurements of the tttt and 
KK decays, allows us to examine this conjecture only 
for /o(1710) which is the only resonance for which we 


have successfully measured branching fractions for both 
decays. We obtain: 


B{fo{171Q)^KK) 
i3(/o(1710) -)> tttt) 


117.6(98) 

37.2(40) 


3.2(4) 


(4) 


The corres pon ding BES result for /(A1710) is i?(BES) = 
2.8(9) [l^ [2^. As quoted in Ref. Q, the Mark III re¬ 
analysis result is i?(Mark III) = 3.7(2®). 

Our result for the ratio R = 3.2(4) differs from 1.33 
expected for a pure glueball by 4.3cr, and can be consid¬ 
ered as clear evidence that /o(1710) is not a pure glue¬ 
ball, although it may well contain appreciable mixture 
of a scalar glueball in its wave function, as many mixing 
models suggest (see also M)- 

To summarize, our measurements have confirmed the 
supernumary nature of light quark scalars, and we find 
that /o(1710) is not a “pure” scalar glueball. Further, 
we find no evidence for the tensor glueball candidate 
C(2230). 

We wish to thank Professors E. Klempt and W. Ochs 
for valuable comments. This investigation was done using 
CLEO data, and as members of the former CLEO Col¬ 
laboration we thank it for this privilege. This research 
was supported by the U.S. Department of Energy. 
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